February 1, 2008 

DIRAC-HARMONIC MAPS 



QUN CHEN, JURGEN JOST, JIAYU LI, GUOFANG WANG 

Abstract. We introduce a functional that couples the nonlinear sigma model 
with a spinor field: L = f M [\d<f>\ 2 + In two dimensions, it is confor- 

mally invariant. The critical points of this functional are called Dirac- harmonic 
maps. We study some geometric and analytic aspects of such maps, in particular 
a removable singularity theorem. 



1. Introduction 

This paper introduces and studies an extension of an established mathematical 
subject, namely harmonic maps from Riemann surfaces into Riemannian manifolds, 
that is motivated by a model from quantum field theory, the supersymmetric non- 
linear sigma model. What distinguishes harmonic maps from surfaces from those 
from higher dimensional domains is the feature of conformal invariance. On one 
hand, this is a global aspect, and it implies that when the domain is the 2-sphere, 
any such map is automatically conformal; thus, the second order equations that 
characterize a harmonic map reduce to first order equations. As is typical in the 
geometric calculus of variations, this leads to important invariants; here, we obtain 
minimal surfaces in Riemannian manifolds as well as quantum cohomology, and 
the theory of pseudo-holomorphic curves in symplectic geometry shares the same 
root. On the other hand, this leads to a non-compact local invariance group and 
turns the regularity and existence problem into a borderline case of the Palais- 
Smale condition. In fact, this conformal invariance connects the local regularity 
theory with global solutions defined on the 2-sphere as we know from the semi- 
nal paper of Sacks- Uhlenbeck [14J. The approach of Sacks-Uhlenbeck depended on 
perturbing the functional to ones that satisfy the Palais-Smale condition and then 
pass to a limit. This, together with a sophisticated local analysis, then yielded 
global existence results, with possible obstructions coming from the second homo- 
topy group of the target. (In the absence of those obstructions, the existence had 
been shown independently by Lemaire Subsequently, alternative existence 

schemes were developed by Struwe[To] and Chang [3] (heat flow method) and Jost 
jH] (local iteration technique). A crucial point of the work of Sacks-Uhlenbeck was 
the removability of isolated singularities through a detailed blow-up analysis. In the 
variational context, it is not difficult to reduce the general regularity question to 
the one of isolated singularities, and therefore, the analysis of Sacks-Uhlenbeck was 
sufficient for the variational existence scheme. The more general regularity question 
for weak solutions of the harmonic map equation was solved later by Helein jjj. 

Now, in the physics literature, the same model goes under the name non-linear 
sigma model; when the target is an (N — l)-dimensional sphere, a case of special 
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interest for quantum field theory, it is called more precisely the non-linear O(N) 
sigma model. Here, our map <f> becomes a Bosonic scalar field satisfying a non- 
linear constraint. Now, this model admits a supersymmetric extension (see j3] for 
a detailed exposition 1 ) where (f> is coupled to a Fermionic field if>. That Fermionic 
field, and then also the Bosonic field, is Grassmann valued, a property needed for 
an additional symmetry exchanging the two types of fields, called super symmetry. 
One can then try to solve the corresponding equations by expansions in the degree 
of the Grassmann variables. 

Here, we rather let the field if) also be real-valued; we then cannot use such ex- 
pansions anymore, but rather get a coupled set of field equations for the two fields 
<f> and Mathematically, <f> as before is a mapping from a Riemann surface into 
a Riemannian manifold, whereas if} is a spinor field with values in the pull-back of 
the tangent bundle of the target under cf>. Thus, with respect to coordinate trans- 
formations on the target, if) transforms as a tangent vector. The important point 
is that the Lagrangian is set up in such a manner as to still be conformally in- 
variant. We thus obtain an extension of the harmonic map model within the class 
of conformally invariant variational problems with a rich internal structure. Here, 
we study the geometric setting of the problem and present a complete variational 
analysis, that is, we prove the removability of isolated singularities and the regu- 
larity of solutions arising from a variational scheme. Our analysis is based on the 
scheme of Sacks-Uhlenbeck |T3] , but because of the coupling between the two fields, 
new difficulties and subtleties arise that make the analysis considerably harder. A 
regularity result in the sense of Helein |7j will be presented elsewhere Having 
thus laid the analytical foundation, the existence question and the construction of 
geometric invariants from the solution spaces can be addressed in subsequent work. 
We should remark that the analysis in our framework is more difficult than in the 
supersymmetric one because we cannot use the expansions of the physicists. On the 
other hand, our present model seems to be very natural from the point of view of 
Riemannian geometry, and it thus falls in a similar category as the Seiberg-Witten 
or Chern-Simons-Higgs models, and this is our main motivation for its study. 

Let us now describe the mathematical framework in more detail: Let (M, h) be 
an oriented, compact Riemannian surface and Pso(2) ~^ S its oriented orthonormal 
frame bundle. A Spin-structure is a lift of the structure group SO (2) to Spin(2), 
i.e. there exists a principal Spin-bundle Ps P in(2) —>■ M such that there is a bundle 
map 

Ps P in(2) ► PsO(2) 

i I 

M — ► M. 

Let E + M := Ps P in(2) x p C be a complex line bundle over M associated to Ps P in(2) 
and to the standard representation p : S 1 — > U(l). This is the bundle of positive 
half-spinors. Its complex conjugate E~M := E+M is called the bundle of negative 
half-spinors. The spinor bundle is EM := S + M © S~M. 



For a variant of this model, see |J 
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There exists a Clifford multiplication 

7 .V x : E+M -> E~M 
7 .V x : X .1/ -> E+M 

denoted by v <S> ip — > i> • V', which satisfies the Clifford relations 

v-w-ip + w- v- ip — —2h(v, w)ip, 

for all i>, iu G TM and ip G SM. 

On the spinor bundle EM there is a hermitian metric (•, •). Let V be the Levi- 
Civita connection on M with respect to h. There is a connection (also denoted by 
V) on SM compatible with the hermitian metric. 

Let be a smooth map from M to another Riemannian manifold (TV, g) of di- 
mension n > 2. Denote <p~ 1 TN the pull-back bundle of TiV by and consider the 
twisted bundle EM <g> 0- 1 TiV . On EM <g> (jr^TN there is a metric induced from the 
metrics on EM and <p~ 1 TN. Also we have a natural connection V on EM®0 _1 TiV 
induced from those on EM and 0~ 1 TiV. In local coordinates, the section -0 of 
EM ® (j)~ l TN can be expressed by 

where is a spinor and {^j} is the natural local basis. V can be expressed by 

n a n a 

V^ = ^V^(x) — (0(a;))+ J2 ^ jk W(x)^(x) — (0(x)). 

i=l y i,j,k=l y 

It is easy to check that 

v(ipi,ip2) = (V^i,^) + (ipuVvfo), 

for any vector field v. 

Now we define the Dirac operator along the map by 

(i.i) pp = + E r} fc v e „^(x)e Q -^(x)^(0(x)), 

i y i,j,k=l y 

where ei, e 2 is the local orthonormal basis of M and ^ := Y^ a =i e <* ' * s the USU& 1 
Dirac operator. The Dirac operator pis formally self-adjoint, i.e., 

(1-2) / ty,#)= / (#,0, 

for all ijj, £ G r(EM <g> (jr^TN), the space of smooth section of EM <g> 0" 1 TiV . Set 

# := {(0, V) | G C°°(M, TV) and G T(EM <g> <fr l TN)}. 
On A", we consider the following functional 

L(0,V) = / [|rf0| 2 + (^,W]v^rf 2 a; 
(1.3) = / [ gi M)h a ^^ + giAm\m j )\^hd 2 x, 

JM ax et OX/3 
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(1.5) 



where r(0) is the tension field of the map (f) and TZ(<ft, ip) G T(4> l TM) defined by 



Here R™j are the components of the Riemannian curvature tensor of h. Note that 
the product (•, •) here is the one of EM. Solutions (0, ip) to ()1.4|) and f)1.5|) are called 
Dirac- harmonic maps. 

As already mentioned, the functional L arises from our study of the supersymmet- 
ric nonlinear sigma model, the difference being that here ip is an ordinary (twisted) 
spinor. Equations (jl.4|) and (jl.5|) couple the harmonic equation and the Dirac equa- 
tion in a rather natural way. 

In this paper, our main aim is to establish some geometric and analytic aspects 
of solutions to equations (jl.4j) and (jl.5|) . We first establish some basic properties 
of the Dirac-harmonic maps and give some examples of nontrivial solutions. We 
also derive some geometric properties of the Dirac-harmonic maps, the conformal 
invariance and the existence of a generalized Hopf differential. Then we prove the 
removablity of singularities for the solutions. 

This paper is organized as follows: in Section 2, we deduce the Euler-Lagrange 
equations, and construct nontrivial solutions; in Section 3, we define the energy- 
momentum tensor of the action L and construct a holomorphic differential (Propo- 
sition 3.2) which plays a role in proving the removable singularity theorem; we also 
establish the basic Weitzenbock formula (Proposition 3.4) for spinor fields satisfying 
(jl.5|) ; in Section 4 we prove the main result about removable singularities (Theorem 



In this section, we establish some basic facts for the functional L and equations 



where r(<p) is the tension field of the map <ft o,nd TZ is defined by M.b\) . 

Proof. Equation (J2.2|) is easy to derive. Consider a family of ip t with di/j t /dt = r\ at 
t = and fix cf). Since pis formally self-adjoint, we have 



(1.6) 




4.7). 



2. Dirac-harmonic maps 



Proposition 2.1. The Euler-Lagrange equations for L are 

(2.1) r(0) = n((f>^) 

(2.2) pip = 0, 




Hence, we get ()2.2|) . 
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Next, we consider a variation {<pt} of </> such that d<fr t /dt = £ at t = and fix ip. 
We choose {e a } as a local orthonormal basis on M such that [e a , d t ] = 0, V ea ep = 
at a considered point. 



(2.3) 



dL{<j> t ) 



U=o 



It is easy to check that 



/ hd<f>t\ 2 \t=o+ [ ^<V,#>|t=o :=/ + //. 



(2.4) J = -2 / r J (0)^ m r- 



See the proof for instance in ^Hj- Now we compute II. First we compute the variation 
of pip. We have 

j t piP = e a -V at V e > 

= e a ■ VeJ? <g> V dt d yt + e a ■ ip* g> V^V^^ 

= e a ■ V e M ® V 8t d yi ) + e« • ^ ® R N (d4>(d t ), d<fr(e a ))d yi . 

Hence, we have 

ii = [ f (ip,j t w)\t=o 



<</>, m ® VaA))U + e« • V* ® R N (d<f>(dt), d<f>(e a ))d yi ) | t=0 



M </M 

TV/ 

M 

V* ® V % a yj )| t=0 + (i/>, e a ■ # ® R N (d(f)(d t ), d<P(e a ))d m )\ t=0 

M 

N, 

M 



(V, e Q • # ® R"(d<t>(d t ), d(/>(e a ))d yi )\ t= o 
{^e a -^ ®R N {rd ym ^ l a d yi )d yi ) 

where we have used (|2.2j) . Consequently, we have 

| t=0 = / [-2^r l (0) + J R m ^(^,V^-^')]f 



M 



and hence (j2~Tjl . □ 

It is obvious that there are two types of trivial solutions. One is (0, 0), where <fi is 
a harmonic map, and another is (y, ip), where y is a point in iV viewed as a constant 
map from M — > N, and ip is a harmonic spinor, i.e., (flip = 0. 

Let us give a construction of non-trivial solutions. Let M = § 2 and N = § 2 . Let 
£§ 2 be the spinor bundle over S 2 with respect to the unique spin structure. For any 
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map : § 2 — > S 2 and a spinor \l/ G r(£8 2 ), we define a spinor field ip along the map 

by 

(2.5) W>,* = e a • ^0*(e Q ), 

where e Q (a = 1,2) is a local basis of T§ 2 . It is clear that Vty,* is well-defined. A 
spinor \1/ G T(S§ 2 ) is called a twistor spinor if 

for any vector field v G T(T§ 2 ). 
Now we have 

Proposition 2.2. Let Vty,* be defined by from a nonconstant map : S 2 — > § 2 

and a spinor Then (0, ip^^) is a Dirac harmonic map if and only if is a 
(possibly branched) conformal map and ^ is a twistor spinor. 

Proof. We first assume that (0, V^,*) is a Dirac-harmonic map. For the spinor field 
if) constructed from (|2.5|) . we always have 

(V^V^'V) = (V0 fc • V0 j ■ V0* ■ *> 
= 0, 

since eg • e 7 = — e 7 • eg for /3 7^ 7 and 0/?0 7 is symmetric. Consequently, 

Hence, if (0, ^,,$) is a Dirac harmonic map, then is harmonic from equation (jl.4j) . 
Choosing suitable coordinates such that V ea e/3 = and x G M, by (|1.5j) at x we 
have 

= 

= e/j ■ V e/3 (e a • ^0*(e a )) 

= e /3 -e a -{V e ^0*(e a ) + W e ^0*(e a )} 

= -{V eQ *0*(e Q ) + v|>T(0)} 

+ei ■ e 2 • {V ei ^0*(e 2 ) - V e2 ^0*(ei) + ^[V ei 0*(e 2 ) - V ei 0*(e 2 )]} 
= -V ea ^0,(e Q ) - ex • e 2 • {V e2 ^0*( ei ) - V ei ^0*(e 2 )} 

Since is conformal (and non-constant), the above equation is equivalent to 

(2.6) et • V ei * = e 2 • V e2 *, 

which is equivalent to the condition that ^ is a twistor spinor. 

Conversely, from the above proof, it is easy to see that if is a conformal map 
and \I/ is a twistor spinor, then (0, , 9 i ) *) is a Dirac-harmonic map. □ 

Remark 2.3. The twistor spinors form a linear space of complex dimension 2, 
whereas the conformal maps form a nonlinear space whose dimension depends on 
their degree; for degree 1, we have dimension 3. 
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3. Geometric aspects of Dirac-harmonic maps 
In this section, we will study geometric properties of Dirac-harmonic maps. First 

Lemma 3.1. The functional L(4>,ip) is conformally invariant. Namely, for any 
conformal diffeomorphism f : M — > M, set 

4> = 4> o / and ijj = X~ l l 2 ip o /. 

Then L(<f),ip) = L((f), tp). Here A is the conformal factor of the conformal map f . 

Proof. The proof is standard. In fact, the following terms are invariant under a 
conformal transformation 

\V(j)\ 2 dvol{g), j \ip,pj))dvol(g) and J \ip\ 4 dvol(g). 

Here we only check the conformality of j (ip, pip)dvol(g). Let g = f*g and <j) the 
usual Dirac operator with respect to the new metric g. By the conformality of /, 
we have g = (X*)g. We identify the new and old spin bundles as in |8}. Recall that 
the relation between two Dirac operators (j) and ^ is (cf. jSj) 

(3.1) ^0 = A~*^, 
from which one can deduce that 

(3.2) lftp = \-%pip, 

and hence the conformal invariance of j (ip, pip)dvol(g). □ 

For a two dimensional harmonic map, there is an important holomorphic qua- 
dratic differential, the Hopf differential. For a Dirac-harmonic map, we also have an 
analogue. 

Let (4>, i/j) be a Dirac-harmonic map. On a small domain Q of M, choose z = x+iy 
a local isothermal parameter z = x + iy with g = ds 2 = p\dz\ 2 . Define 

(3.3) T{z)dz 2 = {{\ct> x \ 2 -\<P y \ 2 -2t(<p x ,<p y )) + (^,d x -V d ^)-t^,d x -V dy ^))}dz 2 . 
Here d x = -i- and d v = 

^ ox y ay 

Proposition 3.2. The quadratic differential Tdz 2 is holomorphic. 

One can prove this proposition by a direct computation, or as a consequence of a 
conservation law as follows. 
Define a two-tensor by 

(3.4) T Q/3 := (2((j) a , (j)p) - <5 a/3 (0 7 , 7 )) + (V>, e a ■ V ep ip), 

where <p a := (f)*(e a ). Here, {e a } is a local orthonormal basis on M and {0 a } a coframe 
dual to {e a }. The tensor T a p6 a ® 13 is called the energy-momentum tensor. Using 
the equation pip = one can verify that T is symmetric: 

T a f3 = Tp a . 
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Proposition 3.3. Let (<f),ip) be a smooth solution of \l-4\j hl.tyl on M, then the 
energy-momentum tensor is conserved, i.e., 

(3.5) ]Tv ea T a/3 = 0. 

a 

Proof. 

V ea T a/3 = V eQ (2(0 a ,0 /3 )-5 a/3 (0 7 ,0 7 )) + V ea (^,e a - V e/3 V) 
:= I + 11. 

We choose a local orthonormal basis {e a } on M such that V ea eg = at a considered 
point. We compute 

/ = 2(V eQ 0*( e a) ) 0*( e /3)) +2(0*( e «)) V 6a 0*( e /3)) 

-25 a/3 (0*(e 7 ), V ea 0*( e 7)) 
= 2(r(0), 4>ts) + 2(0 a , V e/3 0*(e Q )) - 2(0 7 , V e ^(e 7 )> 
= 2(r(0),<M, 

and 

77 = (^ a ,e a • ipp) + (ip,e a ■ V ea V e0 ip) 

Therefore, we have 

(3.6) V e T aP = 2(r(0), <j> p ) + (if,, pfjp). 
Now 

,1 
v 2~ 



2(t(0),^> = 2(-R™W,V<p l -^)d ym ,<P$d yP ) 



= g^R^^-ifS)^ 
(3.7) = RmiijWW 1 -^)<P7- 



We compute pipp = e a ■ V ea V efj if). By a direct computation, we have: 



where Rp M is the curvature operator of the connection V on the spinor bundle EM. 
For this curvature operator, we have (cf. |lUj ) 

(3.8) e a ■ i? SM (e Q , X)^ = -Ric(X) ■ ip\ VX G V(M). 

It follows that 

(if), e a • 7! EM (e a , e^ifj* ® fy) = (tp 4 ® d yj , e a ■ R* M (e a , ® d yi ) 

= g ij (^,e a -R SM (e a ,e p )if; i ) 

= \9iM\ Ric ^a) ■ ^) 
= 0. 
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Therefore 

(ip,flipp) = (1>,e a - V eQ V e/3 ^) 

= (V, V e/3 (e a ■ V e ») + KRjtfpiVtf ■ ^ ® ® ^> 

From (|3.fi)l . (|3.7j) and ()3.9j) we conclude that T Q/ 3 is conserved. □ 

Proof of Proposition 3.2. The proof follows directly from Proposition 3.3. □ 
We now consider (|1.5|) . i.e., 

|ty = + rj. fc (0)0 o ^e Q ■ fc ) ® d yi = 0, 

for : M m -> iV n , ^ e T(SM ® <p~ l TN). Here M n is an n dimensional spin 
manifold with n > 2. Clearly, one can also discuss (1.4)-(1.5) for higher dimensional 
spin manifolds. We call such a ip a harmonic spinor field along the map 0. Note 
that here we do not assume (0, ijj) to be a solution of (jl.4j) . 

For a harmonic spinor field along a map 0, we have the following Weitzenbock 
formula. 

Proposition 3.4. Let M m and N n be Riemannian manifolds, <fi : M — > iV ; and 
= ^ (g g r(EM t&^TN). Then 

(3.9) #V = -V ea V ea ^ + ^ + ^.^(e a • e/J .^*)®A. 
J/ and satisfy hl.ty . then 

(3.10) ^A|^| 2 = |W| 2 + ^/#| 2 - \Rtj M {Vcj> k ■ i>\ V0< ■ i?\ 

where R is the scalar curvature of M , V denotes the connection on EM ® <p~ l TN 
and |0| 2 := ^(0)(^,^). 

Proof. One can apply a general Weitzenbock formula (see for example [TI] or [TUj ) 
to prove the proposition. For the convenience of the reader, we present a proof here. 
Choose an orthonormal basis {e a |a = 1, 2, • • • , m} on M such that V e/3 e a = at a 
considered point. Noting that e 2 = —1, e a ep + epe a = and a ^ ft we have 

(3.11) p 2 i/> = tp ■ e a ■ V e ,V e > 

= -Ve Q Ve Q V' + y^ j e a -e p - R(e a , ep)j) 

a</3 

(3.12) = -V ec , V et > + ^ ^ e a ■ ep ■ R(e a , ep)ip 

a,/3 

where R(-, ■) is the curvature operator on EM <g> <p~ 1 TN, namely, 
Since 

v eQ v e ,0 = (v ea v e/3 # + r} fciP C0^ fc ) ® ^ 
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and 



we have 



R^(e a ,e^ + Rl N <jf a <P^ k ®—. 



^(e a , e/3 )^ = i? SM (e Q , e/3 )#® — + (r;. fciP -r fcpj )C0^ fc ®^- 

Putting this formula into (|3.11|) we have 

p 2 if) = -V ea V ea tP + ^e a -ep-R EM (e a ,e p )ij 

(3-13) ^Ri p ^J p {e a -e p -^)®-^, 

where _R SM (-, •) is the curvature operator on EM. It is known (cf. ^U]) that 

Thus, we obtain 

(3.14) p 2 if) = -V e SeJ + \m + \R kpj Wp{ea-ep-^ k )® A 
from which 



Therefore, 



iA|^| 2 = (V ea V e >,V) + (V e >,V e » 



2 



□ 



Let (N', g') be another Riemannian manifold and / : iV — > N' a smooth map. For 
any (0, if)) E X we set 

(f> = f o <fi and -0' = /*^. 

It is clear that ■0' is a spinor along the map 0'. Let A be the second fundamental 
form of /, i.e., A(X, Y) = (V x df)(Y) for any X, Y G T(TN). The tension fields of 
and 0' have the following relation 

2 

(3.15) r(0') = ^ A(d0(e a ), d0(e a )) + d/(r(0)). 

a=l 

It is also easy to check that the Dirac operators p and p' corresponding to and 
0' respectively are related by the following formula 

(3.16) p'iP' = f,(pp) + A(d<p(e a ), e a • V). 
Here in local coordinates 

5 d 

A(d(f>(e a ), e a -ip) = <p l a A(—, g—)^ • V^, 
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where {^r} is a local basis of N. Furthermore, if / : N — > N' is an isometric 
immersion, then A(-, ■) is the second fundamental form of the submanifold N in N', 
and 

= -P(b X) + Vi^ V X Y = V X Y + A(X, Y) 

VX, Y G r(TiV), £ G T(T- L iV), where P(f ; •) denotes the shape operator. Note that 
(P(f ; X), Y") = (A(X, Y),g). We can rewrite equations (|O j) - (|T3|l in terms of A and 
the geometric data of the ambient space N'. By the equation of Gauss, one has 

(3.17) K{<f>, = P(A(d(f>(e a ), e a ■ </>); V) + ^#(e« ■ V, ^)d<p(e a ). 

Therefore, by using ()3.15|) and ()3.16|) . and identifying with 0' and ■?/> with if}', we 
can rewrite ()1.4|) and (jl.5j) as follows: 

(3.18) r(0) = A(d<t>(e a ), d<P{e a )) + ^R\e a ■ if>, i/>)d<f>{e a ) + P(A(d0(e a ), e a ■ i/j); i/i), 

(3.19) #ty = 4(#(e a ),e a -V). 
In particular, if iV' = R^, then these become 

(3.20) -A0 = A(d0,d0) + P(A(d0(e a ),e Q -V0;V0 

(3.21) ^ = A(#(e a ),e a ^). 



4. Analytic aspects: removable singularities 

Embed (N,h) into some isometrically and denote by A(-,-) the second fun- 
damental form as in the previous section. Any map from (M, g) to (N, h) can be 
seen as a map from (M, g>) to 1R X with 0(x) G X. And any spinor field if) along 
the map can be seen as a X-tuple of (usual) spinors (ip 1 ,^ 2 , ■ ■ ■ ,ip ) satisfying 
the condition that for any normal vector v = 5^ i=1 z^Pj of N at <f>(x), we have 

$>^ = o, 

i 

where {Pj, i — 1, 2, • • • , X} is the standard basis of K^. In this section, we always 
view (0, ip) in this way. By the discussion in the previous section, such a pair (0, ip) 
is a Dirac-harmonic map if and only if (0, ip) satisfies (|3.20)) and (|3.21j) . 

In this section, we will prove the removable singularity theorem for Dirac-harmonic 
maps with "finite energy" defined below. 

Definition 4.1. Let U be a domain on M. The energy of (0,VO on U is: 

(4.0) E{4>,il>,U) := [ (|rf0| 2 + |^| 4 ). 

Ju 

Note that the energy is conformally invariant. It is crucial for our results. 

Before we consider the analytic aspects of Dirac-harmonic maps, let us note that 
on a surface the (usual) Dirac operator ^ can be seen as the Cauchy-Riemann oper- 
ator. Consider M 2 with the Euclidean metric dx 2 + dy 2 . Let e\ = and €2 = ■§- be 
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the standard orthonormal frame. A spinor field is simply a map ^ : R 2 — > A 2 = C 2 , 
and ei and e 2 acting on spinor fields can be identified by multiplication with matrices 



ei 



f 
-1 



e 2 



i 

1 



If * := 



C 2 is a spinor field, then the Dirac operator is 



ftp 



1 
-1 



fdf\ 
\dx/ 



+ 



% 



dg_ 
\dyj 



( — \ 

% 
\ dz) 



where 



d_ 

dz 



d__.d_ 
2 V dx dy 



dx dy 



d__l 

~dl ~ 2 

Therefore, the elliptic estimates developed for (anti-) holomorphic functions can be 
used to study the Dirac equation. 

Proposition 4.2 Let (N n ,h) be a compact Riemannian manifold and (M 2 ,g) a 
surface with a fixed spin structure. Then there is a small constant e > such that 
if (<j), t/>) is a smooth Dirac-harmonic map satisfying 



(4.1) 



(|#| 2 + H 4 )< £ , 



M 



then cf) is constant and consequently if) is a usual harmonic spinor. 

Proof. In view of (13.20}) . we have 

|A0|<C(|d0| 2 + |#|H 2 ), 
where C > is a constant depending only on N. Hence we have 



(4.2) 



HA^y < c(iiv0ii 2 l + 

< C(||V0|| L2 + 

< eC||V0|| £ll | 



i 4 )l|V0|| Ll ,| 



Therefore, if e is small enough, we can show that <fi = const, and hence ip is a 
harmonic spinor. □ 
Now we consider the local behavior of Dirac-harmonic maps. Since they are 
conformally invariant, in the sequel we may assume M to be the unit disk D with 
trivial spin structure. 

Theorem 4.3. There is a small constant e > such that if (0, ip) is a Dirac 
harmonic map satisfying 



(4.3) 

then 
(4.4) 



Jd 



+ 



chd 1 ) + UWchdx) < C{\\V(f>\\ L 2 (D) + W\\n (D) ), 



where C > is a constant depending only on k and the geometry of N. 
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Proof. In the sequel, we denote || ■ Hx^rm by | ■ \o,k, P - If there is no confusion, we 
may drop the subscript D. In the proof, C is a constant, varying from line to line 
and 

Di C D 2 C D 1 C D. 

2 

We devide the proof into several steps. 
Step 1. There is an e > such that 

(4.5) |#bi,o,4 < C{D x ){\d4>\w + IVIm), VD 1 C A 

where C{D 1 ) > is a constant depending only on D 1 . 

Choose a cut-off function 77 : < r] < 1, with 77] pi = 1 and Suppr/ C D. By f!3.20|) 
we have 

|A(#)| < C(|0| + |#|) + |A| oo |#|(| t %0)| + |^|) + 1^1 
< lAUI^H^)! + C(|0| + |#|) + |r/a|, 

where we denote a := P(A(d0(e a ), e a • ip); ip). Thus, for any p > 1, 

(4.6) |A(#)|o, p < |A| oo (|d0||d(#)|) o ,p + C|0|i iP + Mo*- 

Let p = |, and without loss of generality we assume J D (ft = so that |0|i )P < 
C'\d(f)\Q t p, then 

|^|oo(M0|M(#)|) O ,| < |^|oo|#ImI#|o,2, 
from this and ()4.6|) we have 

|#| 2 ,| < C(|^|oo|#|l,4|#|o,2 + |#| ,| + M ,f )' 

By the Sobolev inequality, |^0|i,4 < C'\r](p\ 2 ^, so, 

(4.7) (C'- 1 - ClAU#|o,2)|#|i,4 < C(\d<f>\ 04 + \ V a\ o4 ). 
Moreover, 

IHo,| < C N(\ip\ 2 \vd<t>\)o,l 

= CV(|Vf l#70)-«) Oi| 

< C|^|^|#|l )4 + C^|^4- 

Putting this into ()4.7j) and choosing e small, we get 

|#|l,4 < C(|#| 0j | + V^l#|l,4 + |^lo, 4 ), 

which yields 

Step 2. If e > is small enough, then 

(4.8) Md 2 ,o,,<C(D 2 )H ao ,4, Vg>l, D 2 C D 1 C D, 

where C(D 2 ) > is a constant depending only on .D 2 . 

Choose a cut-off function 77 : < 77 < 1, with 77 1 D 2 = 1 and SuppT? C D 1 . Let 
£ = 777/; have compact support in D l . By the well-known Lichnerowitz' formula, we 
have 
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because the scalar curvature R = on D. Integrating this yields 

I V£ | £,l j0)2 = 1^1^1,0,2 

= \Kv^)\m tQ ,2 = I VV">I> + vW>\dim 

< ^1^1^,0,2 + C\d^\ 2 Mhflfl) 

< C(|Vbi A2 + \d(p\ 2 D i A M 2 D\o,i) 

< C\4j\ D i t0A {i + \d<p\ D i fiA ) 

< C'\iP\ D i m . 

Hence, Step 2 is proved. 

Step 3. If £ is small enough, we have 

(4.9) < C{D 2 )\d(j)\ DA2 , \JD 2 CD 1 CD, 

where C(D 2 ) > is a constant depending only on D 2 . This follows from Steps 1 

and 2. For higher order estimates, it is rather standard. See for example j3j. □ 

From this theorem we know that a sequence of Dirac-harmonic maps with small 
energy has a convergent subsequence. However, if the energy is large, then a blow- 
up may occur. In this case, the concentration of energy may happen. Namely, 
\V(fik(xk)\ 2 — ► oo as k — > oo, for a sequence Xk — > xq as k — > oo. After a suitable 
rescaling, we get a "bubble", an entire solution of (j3.20j) - ()3.21j) with finite energy. 
By the conformal invariance of the Dirac-harmonic map, such an entire solution 
can be viewed as a Dirac-harmonic map from § 2 \{p} — > iV with finite energy. In 
this section, we prove that such a singularity can be removed as in many conformal 
problems. Hence, at the end we obtain a Dirac-harmonic map from S 2 — > N. Using 
the Theorem 4.3, we can describe the behavior of solutions (0, ijj) near a singular 
point as follows: 

Corollary 4.4. There is an e > small enough such that if (0, ip) is a smooth 
solution of jl.^| )-( f7~3j) on D \ {0} with energy E(<f>, ip, D) < e, then for any x G Di 

2 )l 



(4.10) \d<f>(x)\\x\ < C( \d(j)\ 

' D(2\x\) 



(4.11) \i)(x)\\x\* + \Vip(x)\\x\3 <C( I I VI 

D(2\x\) 



I 1^ 



4 



Proof. Fix any xq G D \ {0}, define <fi and ip by 

4>{x) := 4>(x + \x \x) and ^(x) := \x \^ip(x + |x |x). 

It is clear that (0,V>) is a C°° solution of (jl.4j) - (fT3j) on D and E((j),ip,D) < e. 
Applying Theorem 4.3, we have 

|#|l°°(Di) < C|#b,0,2 

2" 

and 

Mci(£>i) < C|#b,0,4- 
2 
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Scaling back, we prove the Corollary. □ 

Lemma 4.5. Let (4>,ip) be a C°° solution of ji.^ -j TOp on D \ {0} satisfying 
E(<f),i/j,D) < e, then, 

1 1 ' 12, 



r2 ./I! 



(4.12) / -|<^| 2 d# = / \(f> r \ 2 d6+ / (il>,d r -1> r )d9 

) Jo 

»2tt /.2tt i 

V| 2 rf0 — / -(ip,d e -ip e )d0, 
o Jo r 

where (r,9) are the polar coordinates in D centered at 0, ip r := Va r ip. 

Proof. From Proposition 3.2, 

T = [|0,| 2 -|0j 2 -2z(^,^)] + [(^,^-V 9 ^)-^(V',^-Va^))] 
=: (A-iB) + (A'-iB'), z = x + iyeD 

is holomorphic on D\ {0}. By Corollary 4.4, we know that 

\A-iB\<2(\<p x \ 2 + \<P v \ 2 )<C\z\- 2 . 

Noting that V ' e j\) = (V^* + T)^^) <g> d yi , we have 

|W| <C(|W| + |d0|H)- 
By Corollary 4.4 again, we have 

\A' - iB'\ < 2M|VV>| < CdVllVVI + |<ty||Vf) < C\z\~ 2 . 
Therefore |T(z)| < C\z\~ 2 . Furthermore, 

/ \A'-iB'\ < 2 / IVHWI < oo 

and / |A — iB\ < 2 J D \d(p\ 2 < oo. Thus, j D \T(z)\ < oo, which implies that T(z) 
has a pole at 2 = of order at most one. Hence, zT(z) is holomorphic on D and 

r /*27r 

(4.13) = lm[/ zT(z)dz] = Re[z 2 T(z)}d9. 

J\z\=r Jo 

It is easy to compute that 

(4.14) Re[z 2 T(^)] = r 2 [(A cos 26 + B sin 26) + (A' cos 26> + B' sin 20)] 

= r 2 |0 r | 2 -|0,| 2 -(V^-^>. 
Now the desired equalities follow. □ 

Remark 4.6. Integrating (|4.12|) yields: 

(4.15) [ \<p r \ 2 - [ \\M 2 = ~ I MA- 



D JD r2 JD 



From (I4TT21 and 



/ I0,| 2 + / \\U 2 = [ m 2 := E^) 

J \z\=r J \z\=r J\z\=r 
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we have 

(4.16) / \(f)r\ 2 = \e t + h 



2" r ' 



and 



1 ■ ■ ,2 1 



r ■ 



(4.17) = _*-_,, 

where I r := - f\ z \ =r (i>, d r ■ ip r ). 
Now we can state the following 

Theorem 4.7 (Removable singularity theorem) Let ((f), ip) be a solution of ( |1.^| ) 
and M.5)) which is C°° on U\ {p} for some p G U C M. If (0, ip) has finite energy, 
then (<f>,if>) extends to a C°° solution on U . 

Proof. By rescaling, we may assume that 

(4.18) / (|d0| 2 + H 4 )<e. 

JD{2) 

Choose a function q(r) on D which is piecewise linear in logr with 

g(2- m ) = ^- 1^ <p(2- m ,e)de, 

then we have (cf. [T4*] ) 

(4.19) / \dq-d<P\ 2 = [ (q-<j>)<j> r - [ (q ~ <f>)A(q - <f>) 

JD Jr=l JD 

with 

(4.20) \q - 0U := \q - <j>\ c o {D) < 2 3 v^. 
Using (3.21), we have A(q — (f>) = —A(f> = —A(<f>)(d<f), d<p) — a, so, 

(g-0)A(g-0)| < |g-0|oo|^U / \d(j)\ 2 + \q - 0U / M 

D JD JD 

(4.21) < 5 [ \dcj)\ 2 + Cy/i [ \^\ 2 \d(j)l 

JD JD 

where 2 3 y / £| J 4| 00 < 5 for some constant 5 > small, and C > is a constant. As 
for the first term on the RHS of (jOTSjl . 

(9-0)0, < (/ \q-<t>\ 2 )hf I0,l¥ 

r=l J r=l J r=l 



(4.22) < 
On the other hand, by (4.16), 



r=l J r=l 



(4.23) / \dq-d<P\*> / = -J, 



D JD 
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where I := - f D (ip, d r ■ ip r ). Inserting ijOTl) . KT2h and (H^l) into (14*131) . we get 
1 



2 



D * Jr=l Jr=l 



-I < (I \U 2 )H \4>A 2 )> 



(4.24) +5 / \d(f>\ 2 + Cy/i / |V| 2 |#I- 

JD JD 

By (jjHEI and (jQ7jl . we have 



r=l J r=l 



< -E x = - 
- 2 2 



r=l 



It follows 



(4.25) (1-25) / |d0| 2 < / |d0| 2 + [2C\/i / / (V,9r"0r)]. 

Jr=l JD JD 

By a scaling argument, this yields 
(1-25) [ \d(P\ 2 < r[ \d<l>\ 2 + [2Cy/l [ |^| 2 |#|-/ ty,d r "*l>r)\ 

J D r JdD r JD r J D r 

< rf |#| 2 + CVi/ W + Cy/if \d<P\ 2 +f |^||V^| 

(4.26) < r/ \d<f)\ 2 + C^[ \d<p\ 2 + C f \^ + C f |W|i 

Now we need the following lemma. It is the elliptic estimate with boundary. For 
completeness we give a proof here. 

Lemma 4.8. Let u be a complex function satisfying 

u 27) { ® u = f> in D > 

\ u\ dD = <p, 

with ip e L 1,p (dD) and f G L P (D) for some p > 1, where D is the unit disc on M 2 
centered at the origin, then the following estimate holds 

(4.28) W\D,x,p < C(\f\ Dfi , P + \<p\dD,i, P )- 

If instead u satisfies 

du = f, in D, 



(4.29) 

then the same estimate holds. 



U\8D = <P, 



Proof. We first consider the following boundary value problem: 

dwi = 0, 

Rewxldc = (p, J 9D Imwi = 2nc , 



(4.30) 



where cq := j- f dD limp. It is clear that (cf. [2], Theorem 38): 
(4.31) \wi\c*(D) + |Vwi| LP(D ) < C(\(p\ dDilj , + c ) 

where C > is a constant. 
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Next, let u>2 :— u — Wi, then it satisfies 

( 4 - 32 ) | 9w 2 = / in D 

[Rew 2 \dD = 0, J dD lmw 2 = 0. 

The elliptic estimates yield that 

(4.33) | W2 1 £!.*>(£>) < C\f\ LP{D) . 

In fact, by the Schwarz-Poisson formula (see Theorem 21 in |2.), one has 

(4.34) t „ 2 = _I^ i( _Z_ + _iL ?Rdl) . 

From this and the boundedness of the Riesz transformation, we obtain (|4.33|) . Com- 
bining (|4.31|) and ()4.33|) then gives the desired estimate ()4.28|) . The proof is similar 
if u satisfies (PT2H . □ 
Now we return to the proof of Theorem 4.7. Recall (3.22) 

(4.35) #0 = A(d0(e a ),e a -<0) in D \ {0}. 

We choose a cut-off function r] £ e C^(D 2£ ) such that r] e = 1 in D £ (0) and |<ir7 e | < 
C/e. Then we have 

- = (1 - Ve)A(d(j)(e Q ) } e a ■*/))- V Ve ■ 0. 

From Lemma 4.8, we have 

(4.36) |(1 - riM D>lt * < C\d<p\ D ,oMD,oA + C|V>U,i,f + C\V Ve • ^| da| . 
Letting e — > 0, using 

lim-V / = 0, 



the smallness of |g?0|z?,o,2 an d the Sobolev embedding theorem, we obtain 

0l¥<C(/ |V#)i+C(/ H 4 )i 



By rescaling, we have for any < r < 1 



)2 < C(r / |W|3)3 + C(r 

D r J0D r JdD r 



Thus, 

(4.37) f \ip\ A <Cr( \Vip\%+Cr[ |0| 4 . 

JD r JdDr JdDr 

Let 4> '■= ^ Ibd^- Note that 

^-^)=A(#(e a ),e a .(V'-^))+A(d0(e a ),e a -^) in -D \ {0}. 

By an argument similar to the one used in obtaining ()4.36|) and using the Poincare 
inequality, we have 

(4.38) +C\V^\ d n oA . 
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Again, by the smallness of |£Z^»|d,o,2 we obtain 



'D JdD JD 



< c(l |w>|S)* + c(/ 

I'dD JdD JD 



So we have 



|vv>|3 < c |v#r + c(/ |^| 4 )3( 

D JdD JdD JD 

< C I |V#+£! / |rf0| 2 + - / 
JdD JD £ l Jd 



2\i 



M 4 , 

3D 



where S\ > is a small constant. Hence, for < r < 1, 

(4.39) / |V#<Cr/ |VV|5+ei/ \d<P\ 2 + —f |i/>| 4 . 

Putting ()4.26|) . ()4.37|) and ()4.39|) together, we have for any < r < 1 and some 
constant C > 

2 i l" L/.I4 i / lr7„/.ll ^ ri„( I I J J. 12 



(4.40) / \d<p\>+ / |^| 4 +/ |V-0ia <Cr(/ |d^+/ |^| 4 +/ |V^|3). 

JL> r JZ) r J D r J 8D r J3D r JdD r 

Denote F(r) := f \d(j)\ 2 + f |i/>| 4 + / \Vif>\*. (IQOl implies that 

(4.41) F(r) < CrF'(r). 
Integrating this inequality yields 

(4.42) F(r) < F(l)rv. 

From this we can easily conclude that there are some (3 > 1 and 2 > q > | such 
that 

(4.43) E L l 'W{D), i> e L l,q (D). 

Now let us consider tp. Recall that (f)ip — A(d(j)(e a ), e a ■ ip) := a, and note that 



\a\ qi < C / (\d<PM) 



\9l 



D i JDx 

2 3 



< c{[ i#r')^(/ 

JD, Jd 



1919 



7 2 



where q\ > 0,p' and </ are constants chosen as follows 

20q 2q 2q 



9l = 7775 775 TO P 



2^-09-l)g' " 2g-(2- ff )gi' * (2-g)gT 
Since /3 > 1, we have qi > q and 

(4.44) / \a\ qi <C{[ \d(f)\ w y(f < oo, 



JDi JD 1 JD X 

2 2 2 

where q* := We note that ib G L qi because of q\ < q* . 

From the regularity theory of Cauchy-Riemann operators we have 

(4.45) ip 6 Lf with , 1 = 2p _\ P ;_ 1)q >,- 



20 
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From above, we see that <fi G L l,2 @ and ip G L 1,q (4/3 < q < 2) imply G L 1 ' 2 ^ 
and tp G L 1 ' 91 . By iteration, G L 1,2/3 and -0 G L 1,9 " imply G L 1 ' 213 and -0 G L 1,9n+1 
with 

Wl 2/3-(/3-l)^' 

Since 

gn+i _ 2/3 3/3 

g n 2(3-{(3-l)q n (3 + 2 
there exists g„ such that g n > 2. Therefore, we have that <fi G L 1,Po (D) and -0 G 
L 1,q °(D) for some p > 2 and q > 2. We can then conclude that ((f), tp) is smooth 
on D through the standard bootstrap method. We omit the details here. This 
completes the proof of the theorem. □ 
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